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The vertex models have been intensively studied and 
their thermodynamic properties are well known inside 
the solvable parameter area. 1,2 ) On the other hand, out- 
side the solvable area the model is not fully analyzed. 
Takasaki et al. numerically investigated a vertex model, 
which allows 7 vertex configurations and contains 2 pa- 
rameters. 3 ' They observed a phase transition that be- 
longs to the Ising universality class. In this article we 
study a symmetric vertex model, that allows 10 ver- 
tex configurations, by use of the corner transfer matrix 
renormalization group (CTMRG), 4 ) a variant of the den- 
sity matrix renormalization group. 5 ~ 7 ^ As we report in 
the following the model has a critical point. 



where P e is the probability to observe vertices with £ 
numbers of 1, is the increasing function of x. 
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Fig. 1. Allowed vertex configurations and their Boltzmann 
weights. The open and black circles, respectively, denote spins 
that take and 1. The shaded area will be used when we plot 
the spin configuration snapshots in Fig. 2 and 5. 



Let us introduce the vertex model that we study in 
the following. Consider a square lattice, where there is 
a 2-state spin variable (a = or 1) on each bond. Thus 
a lattice point is surrounded by 4 spins. We impose a 
local constraint for these 4 spins as shown in Fig.l. We 
assign the following Boltzmann weights 



W = l, W 1 =x, W 2 =ax 2 , W 4 = x 4 



(1) 



for these 10 local configurations, where the subscripts 
denote the sum of 4 spin variables. The parameter x is 
positive. The expectation value of a spin 



Fig. 2. Spin configuration snapshot at a = and x = 1.23. For 
each vertex, the shaded area in Fig.l is drawn. This configu- 
ration snapshot is obtained by recently developed extension of 
CTMRG. 8 ' 9 ' 



When the parameter a is 0, only 6 vertex configura- 
tions are allowed. Under this strong constraint, rect- 
angular areas inside which all spins are a = 1, appear 
separately in the 'sea' of a = 0, as shown in Fig. 2. The 
expectation value p(x, a = 0) gradually increases with x 
till x = 1.23. At this point the system shows first order 
phase transition, and p jumps to unity. (See Fig. 3.) 




p(x 1 a) = (a) = (P 1 +2P 2 +AP 4 )/4. 



(2) 



Fig. 3. Expectation value p(x, a) for a = 0, 0.3599, and 0.6. 



In contrast, when a = 0.6 there is no discontinuity 
nor singularity in p(x, a = 0.6). Thus there should be a 
critical point at certain values of a and x. We trace the 
line of the first order transition in the parameter space 
as shown in Fig. 4. The critical value of a is roughly 
determined as 0.36. We further perform scaling analysis 



2 



Kouji Ueda, Ryota Otani, Yukinobu NlSHlO, Andrej Gendiar, and Tomotoshi NlSHlNO 



for p(x, a), and find one of the parameter line 

a = -0.2736(a; - x c ) + a c (3) 

shown in the inset, that passes the critical point at 
(x c ,a c ) = (1.31438,0.3599). Figure 5 shows the con- 
figuration snapshot at the criticality. 
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Fig. 4. Phase Diagram. The line of the first order phase tran- 
sition ends at the critical point (x c ,o c ) = (L3f438, 0.3599). 
Cross marks represent data points used for the snapshot obser- 
vation in Fig. 2 and Fig. 5. Inset: the parameter line used for the 
scaling analysis in Fig. 6. 




Fig. 5. Snapshot at the criticality. 



In order to determine the spin expectation value 
p(x c , a c ) at the criticality, we assume a scaling form 



p(x, a) - p(x c ,a c ) oc \x c - x^ 



(4) 



on the parameter line that passes through the critical 
point, 10 -' and plot p(x, a) on the parameter line given by 
Eq.(3) with respect to a various power of \x c — x\. As 
a result we confirm that p(x, a) shows linear dependence 
with \x c — x] 1 / 15 as shown in Fig. 6, where the two fitted 
lines are represented by the following equations 

0.369(3; -a; c ) 1/15 + 0.716 (x > x c ) 

-0A22{x c - x) 1/15 + 0.716 (x < x c ) . (5) 

In conclusion the exponent p. in Eq.(4) is 1/15 and the 
value p(x Cl a c ) is 0.716. Thus the critical point belongs 
to the same universality class as the 2D Ising model, 
whose magnetization at critical temperature obeys the 



scaling form 



where S = 15. 



M(T c ,h) ocH 1 /* 



(6) 
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Fig. 6. Spin expectation value p on the parameter line given by 
Eq.(3). The exponent fj, is fixed to 1/15. 



To summarize, we investigated the phase transition of 
a vertex model with 10 local configurations, and found a 
critical point that belongs to the Ising universality class. 
A possible extension of the current study is toward the 
direction of the lattice polymer. 11,12 ) For example, intro- 
ducing two new vertices to those vertices shown in Fig. 1 , 
one obtains a unified model with the 7-vertex case stud- 
ied by Takasaki et al. that corresponds to a straight line 
polymer in two dimension. 3 ) 
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